High- Frequency Properties of a Graphene Nanoribbon Field-Effect Transistor 
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We propose an analytical device model for a graphene nanoribbon field-effect transistor (GNR- 
FET). The GNR-FET under consideration is based on a heterostructure which consists of an array 
of nanoribbons clad between the highly conducting substrate (the back gate) and the top gate 
controlling the dc and ac source-drain currents. Using the model developed, we derive explicit ana- 
lytical formulas for the GNR-FET transconductance as a function of the signal frequency, collision 
frequency of electrons, and the top gate length. The transition from the ballistic and to strongly 
coUisional electron transport is considered. 
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PACS numbers: 73.50.Pz, 73.63.-b, 81.05.Uw 



I. INTRODUCTION 

The utilization of the patterned graphene which con- 
stitutes an array of sufficiently narrow graphene nanorib- 
bons as the basis for field-effect transistors (GNR-FETs) 
opens up wide opportunities for the energy band engi- 
neering promoting the achievement of the optimal de- 
vice parameters. By properly choosing the width of the 
nanoribbons, one can fabricate the graphene structures 
with relatively wide band gap but rather high electron 
(hole) mobiHty [l|,[2 (see also Refs. [ii,[i,iJ3 and the 
references therein). 

In this paper, we present an analytical device model 
for a graphene-nanoribbon FET (GNR-FET) and obtain 
the device characteristics. The GNR-FET under con- 
sideration is based on a patterned graphene layer which 
constitutes a dense array of parallel nanoribbons of width 
d with the spacing between the nanoribbons ds ^ d. The 
nanoribbon edges are connected to the conducting pads 
serving as the transistor source and drain. A highly con- 
ducting substrate plays the role of the back gate, whereas 
the top gate serves to control the source-drain current. 
The device structure is schematically shown in Fig. 1. 
Using the developed model, we calculate the potential 
distributions in the GNR-FET as a function of the back 
gate, top gate, and drain voltages, VJ,, Vg, and Vd (reck- 
oned from the potential of the source contact), respec- 
tively, and the GNR-FET dc characteristics. This cor- 
responds to a GNR-FET in the common-source circuit. 
The case of common drain will briefly be discussed as 
well. For the sake of definiteness, the back gate voltage 
Vb and the top gate voltage Vg are assumed to be, re- 
spectively, positive and negative [Vb > and V^ < 0) 
with respect to the potential of the source contact, so we 
consider GNR-FETs with the channel of n-type. 

The paper is organized as follows. In Sec. II, the GNR- 
FET device model is considered. Here the main equa- 



tions governing the GNR-FET operation are presented. 
Section III deals with the derivation of the GNR-FET 
steady-state characteristics. The results of this section 
generalize those obtained previously Q for the case of 
ballistic electron transport for arbitrary values of the 
electron collision frequency. In Sec. IV, we derive a gen- 
eral formula for the ac source-drain current as a function 
of the signal frequency and the device material and struc- 
tural parameters. Section V is devoted to the analysis 
of the GNR-FET transconductance in a wide range of 
the signal frequencies and device parameters. Section VI 
deals with the interpretation of the obtained character- 
istics and discussion of their limitations. In Sec. VII, we 
draw the main results. 
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FIG. 1: Schematic view of (a) a GNR-FET structure and (b) 
its band diagram at the applied voltages. 



II. GNR-FET MODEL 



Here 



We consider a GNR-FET with the structure shown in 
Fig. 1(a). The back gate and drain voltages T4 and V^ 
(with Vd < Vg) are positive. As a result, the electron 
densities in the channel sections adjacent to the sourse 
and drain (not covered by the top gate) are, respectively, 
given by [7| 
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(1) 



where ae is the dielectric constant, Wt is the thickness of 
the layer between the GNR array and the back gate. The 
top gate voltage comprises the dc and ac components, Vg 
and SV{t), respectively (Vg < 0, \SV{t)\ < |Vg|). When 
the absolute value of the dc component of the top gate 
voltage Vg is sufficiently large {\Vg\ > VbWg/Wb), the 
channel section beneath the top gate is essentially de- 
pleted and potential barriers for the electrons incident 
from the source and drain sections of the channel are 
formed (see Fig. 1(b)). In the case of a GNR-FET with 
sufhciently long top gate {Lg :^ A = y/WbWgJS), the 
height of the barriers can be presented as Q 
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(2) 



where Aq = 

A'* 



eiVgWb + VbWg)/iWb + Wg) nud Ag = 
eVd- As seen from Eq. (2), the variation of the top 
gate voltage leads to the variarion of the potential barrier 
heights and, consequetly, to the variation of the electron 
current above the barrier. 

We assume that the electron energy spectrum of the 
lowest GNR subband can be presented as 
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dp 
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is the electron velocity with the momentum p, Vp is the 
collision frequency of electrons associated with the disor- 
der (including edge roughnesses) and acoustic phonons. 
We disregard the inelasticity of the scattering on acous- 
tic phonons, so that all the scattering mechanisms under 
consideration result in the change of the electron mo- 
mentum from p to —p. The probability of the electron 
elastic scattering is essentially determined by the density 
of state. Taking into account coUisional broadening [^, Q 
(see also, Refs. [iOilllI; we approximate the dependence 
in question as Vp = v\/"f'^ + 1/1/7^ + {p/ptY, where v 
is the value of the coUisional frequency of thermal elec- 
trons, 7 chracterizes the collision broadening, and pT = 
\/2kBTm*. If \p\/pT < 7, fp — v\/l'^ + 1/7 — const, 
whereas at \p\/pT ^ 1 (but \p\ < p^ 
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The bondary conditions for f{p,x,t) are set at the 
points in the electron channel beneath the edges of the 
top gate, i.e., at a; = and x — Lg: 
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Here 
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Here the energy is reconed from the subband bottom, 
A is the energy gap, p is the momentum in along the 
nanoribbon, m* = A/2w^ is the electron effective mass 
(near the subband bottom), and vw — 10^ cm/s is the 
characteristic velocity of electrons in graphene (and in 
GNRs). The energy gap A and the effective mass m* 
depend on the GNR widths d. In the simplest model, 
A c>c m* (X d^^. The electron transport in the section of 
the electron channel beneath the top gate (the gated sec- 
tion), in which the electron density is small, so that one 
can neglect the electron-electron scattering, can be de- 
scribed by the following kinetic equation for the electron 
distribution function: 
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are the distribution functions of the electrons incoming 
to the channel section in question from the source and 
drain sides, respectively, ks is the Boltzmann constant, 
and T is the temperature. 



III. STEADY-STATE CHARACTERISTICS 



Solwing Eq. (4) with boundary conditions (6) at 
SVg{t) — 0, for the dc component of the electron dis- 
tribution function fa{p,x) we obtain the following rela- 
tionship: 
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Here /3p — VpLg/vp. For definiteness, here and in the fol- 
lowing we assume that the modulus of the top gate volt- 
age |V^| is sufficiently large, so that the channel beneath 
this gate is essentially depleted Q ■ At such top gate volt- 
ages, the GNR-FET dc transconductance is much larger 



that when iVgj is small. 



The source-drain dc current can be calculated using 
the following formula: 



Jo(x) = e / dpvp [fo{p, x) - /o(-p, x)] ^ Jq ^ const. 
Jo 

(8) 
Considering Eqs. (7) and (8), we obtain 
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As follows from Eq. (3), 
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Since the kinetic energy of electrons propagating in the 
gated section of the channel of realistic GNR-FETs can 
be assumed small in comparison with the energy gap 
(Ej, ~ fcsT <C A), from Eq. (11) one obtains the fol- 
lowing simple relation: v{ep) ~ ^2ep/m*. Taking this 
into account, from Eq. (9) we obtain 
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Here 
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is the dc source-drain current in the ballistic regime [7|, 
and the factor 



Go = 
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descri bes the ef fect of electron collisions. Here /3 = 
V Lg^Jm* l2kBT can be called the ballistic parameter 
This para meter can also be presented as /3 = vt, where 
T — Lg^m* jlksT is the characteristic electron tran- 
sit time beneath the top gate and VT\Jnn* jlksT is the 
thermal electron velocity. The case /3 <C 1 corresponds to 
near ballistic transport of the majority of electrons (ex- 
cept, possibly, a fraction of low energy electrons). In the 
opposite case (/3 ^ 1), the collisions substantially affect 
the electron transport. When /3 tends to zero, Gq turns 
to unity. Hence, in the ballistic regime of the electron 
transport through the section of the channel covered by 
the top gate, Eq. (11) yields Jo = J(f , where the value 
J(f as a function of the gate and drain voltages coincides 
with that obtained previously 7] . In the case of collision- 
dominated electron transport (/3 3> 1), from Eq. (13) we 
obtain 



Gn '^ — 



(14) 



where a is a coefficient (weakly dependent on 7). 

Equations (11) - (14) desribe the GNR-FET dc 
current- voltage characteristics, i.e., the dependences of 
the dc source-drain current on the gate voltages as well 
as the drain voltage, which generalize those derived pre- 
viously 7] for the purely ballistic transport. 




FIG. 2: Frequency dependences of the GNR-FET transcon- 
ductance modulus for fiflerent values of ballistic parameter /3. 
Inset shows Go as a funcion of this parameter. 




where S V^ and oj are the signal anipHtudc and frequency. 
In such a case, Eq. (4) for the ac component of the distri- 
bution function 5 /^ {ja, x) in the channel section beneath 
the top gate reduces to the following: 

(i O T ( 7} T") 

-~iujSf^{p,x)+Vp ^ — ■ — = ~i^p[6f^{p,x)-Sfu;{-p,x)]. 

(15) 
As follows from Eq. (6), the bondary conditions for 
5 fu]{p,x) in the most practical case eVd S> fcsT can be 
presented in the following form: 



FIG. 3: The GNR-FET amplitude-phase diagram for the 
same parameters as in Fig. 2. 
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IV. SMALL-SIGNAL ANALISIS 



SfLu{P,x)\p<o,x=Lg = 0. 



(16) 



Consider now the case when the top gate voltage com- 
prises a small-signal comoncnt SV{t) — 5V^ exp(— ia;i), 



J 



Solwing Eq. (15) with boundary conditions (16) and 
using Eq. (2), we arrive at the following formula: 
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where Co — yju){u} + i2vp). Taking into account that the ac component of the electron current 5 J^ in the channel 
beneath the top gate is given by 



SJ^{x)^e / dpvp[d f^{p,x) - 6 fuji-p,x)] 
Jo 



and using Eq. (20), we obtain 
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In particular, near the drain-side edge of the top gate {x — Lg), Eq. (19) yields 

f /'OO - 

SJu>iLg) = \ / dpvpfMip)—— — -^ / --T / ^ — F ^2 r~T i \ 

tJo {uo + uoY(syi^{-iu}LglVp\- {uj -LoYcx^{iuoLglVp) 



(19) 



X exp 



eiVbWg + YgWb) 



{Wb + Wg)kBT 



Wb Ae'SV^ 

{Wb + Wg) kBT 



r 



(20) 



Equation (20) can be presented as 

,B Wb e5V^ 



where 
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FIG. 4: Normalized transconductance vs signal frequency 
/ = cj/27r calculated for ballistic (/3 = 0) GNR-FETs with 
different top gate lengths. 



FIG. 6: Normalized threshold frequency utT vs ballistic pa- 
rameter j3. 
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FIG. 5: The same as in Fig. 4 but for GNR-FETs with essen- 
tial electron scattering {13 = 1). 
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. (22) 
where Jq is given by Eq. (12) and Cj = 

In this case, setting ly = in Eq. (22), we obtain 



G^= j d^ expl ?— = -^ 



Naturally, at w = 0, Eq. (23) yields G^ = Gq == 1. 



(23) 



V. TRANSCONDUCTANCE 

The electrons passing from the the channel section be- 
neath the top gate into the drain section of the channel 
between the edge of the top gate {x = Lg) and the edge 
of the highly conducting regions adjucent to the drain 
contact (x = Lg -\- L^), induce the current through the 



latter. This induced source-drain terminal ac current can 
be calculated using the Shockley-Ramo theorem [ij, [lj| 
accounting for the features of the geometry of all highly 
conducting contacts [IJ, |la|- The terminal ac current 
5 J^ is given by 



5 Jul — 6 J^{Lg)D^ — Jq 



Wb 



5<5K. 



{Wb + Wg) ksT 



G^D^. (24) 



Here Ld is the length of the depleted region between the 
top gate and the drain (see Fig. 1) determined by the 
applied dc voltages. In particular, if Vd > Vb, this de- 
pleted region extends to the drain contact. The electron 
drain transit-time factor D^ can be calculated using the 
following formula: 



''-f 
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Lg+Ld 



dx g{x) exp 
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where Td is the electron transit time across the gate-drain 
region, and .9(2;) is the form- factor accounting for the 
shape of the highly conducting regions (top gate and 
drain contact as well as the higly conducting portion 
of the channel near the latter contact). The electric 
field in the gate-drain region, so that their average ve- 
locity can substantially exceed the electron thermal ve- 
locity. Hence, one might expect that the electron tran- 
sit time across the gate-drain region is sufficiently short 
Td <C T. As a result, in the most interesting frequency 
range lot ^ 1, one obtains LOTd <C 1. In this case, one can 
put D^ = 1 . The transit-time effects which lead to a sig- 
nificant frequency dependence of D,^ and, in particular, 
to a marked decrease in | _D„ | at elevated frequencies are 
discussed in the following. 

Figure 2 shows the frequency dependences of the GNR- 
FET transconductance modulus |^a;| = \GuiD^\ ~ \Gui\ 
calculated for different values of the ballistic parameter 
(3. It was assumed that 7 = 0.5. The inset shows Go as 
a function of /3. As it was demonstrated the Go vs (3 de- 
pendence is almost insensitive to variations of parameter 
7; the maximum deviation from the dependence shown 




FIG. 7: Drain transit-time factor vs frequency: (a) for bulk 
drain contact and (b) for blade-like drain contact. 



in the inset in Fig. 2 was less than 3.5% when 7 varied 
from 7 = to 2. In this and the following figures, we 
set A = 0.4 eV and T — 300 K. Figure 3 demonstrates 
the amplitude-phase diagram calculated for the same pa- 
rameters as in Fig. 2. The normalized transconductance 
\Guj\/Go as a function of the signal frequency / = Ci;/27r 
calculated for GNR-FETs with different top gate lengths 
is show in Figs. 4 and 5. 

Figure 6 shows the dependence of the normalized 
theshold frequency uJtT at which \Gui\/Go = l/\/2- It is 
instructive that the obtained dependences exhibit max- 
ima at certain values of the ballistic parameter (3 (not at 
f3 — a.s one might expect). Moreover, the height of these 
maxima decreases with increasing parameter 7, which 
characterizes the coUisional broadening (see Sec. II). 



VI. DISCUSSION 

The effect of the dependence of the threshold frequency 
ujt on parameters /? and 7 can be attributed to the fol- 
lowing. First, at the ballistic electron transport, the fre- 
quency dependence of G^j is actually determined not only 
by the characteristic transit time tau but also by the ther- 
mal spread in the electron velocities [l^, [l3| • The latter 
is about the thermal electron velocity vt- The scatter- 
ing of electrons, particularly strongly dependent on the 
electron momentum, can influence the effective spead in 
velocities. Indeed, when 7^1, i.e., when the low-energy 
electrons (with small |p|) exhibit rather strong scattering, 
their contribution to the current decreases with simul- 
tanious decrease in the effective spead in velocities. This 
can be seen, in particular, from Eq. (7) in which the 
asymmetric part of the distribution function /o {p, x) is 
fo{p, x) - fo{-p, x) ex /m(p)/(1 + Pp). Since the term /3p 
sharply increases when p tends to zero, this function can 
be fairly small in the range of small \p\ (especially at small 
7), while at large |p| it tends to the Maxwellian func- 
tion. This implies that, the momentun distribution can 
shrink and the velocity spread can decrease due to spe- 



cific feature of the electron scattering in GNRs. A signif- 
icant increase in the threshold frequency can be achieved 
if the directional electron velocity markedly exceeds the 
thermal velocity. This can occur in GNR-FETs with the 
GNRs with varied widths (increasing from the source-side 
edge of the top gate to its drain-side) and the pertinent 
decrease in the energy gap. In such hetero-GNRs, the 
electrons can accelerate propagating from the source to 
drain (see, for instance, Ref. 18]). 

As was mentioned above, factor Z3^ associated with 
the electron transit-time effects in the region between 
the top gate and drain contact depends on the tran- 
sit angle ujTd and the form- factor g{x). Figure 7 
shows the frequency dependence of \D^\ for the fol- 
lowing two cases: (a) g{x) = 1 and (b) g{x) = 
Ld^/[x - {Lg - Ld/2)2 - Ly4\ [11, [H. The form-factor 
(a) can be attributed to the case when the drain con- 
tact is rather bulk, so that the source-drain ac current 
is induced primarily in this contact (not in the highly 
conducting portion of the channel adjacent to the drain 
contact). The form- factor (b) corresponds to the case 
when the ac current is induced mostly in the highly con- 
ducting region or when the drain contact is a blade-like. 
Both such cases can take place depending on the spac- 
ing, Lgd, between the top gate and the drain contact 
and the bias voltages Vb and Vd- At Vd > Vb, the re- 
gion of the channel between the top gate and the drain 
contact is virtually depleted. Hence, if the drain con- 
tact is bulk, g(x) ~ 1 with Ld = L^d- The frequency 
dependences of \Di^\ calculated for the cases under con- 
sideration are shown in Fig. 7. As seen from Fig. 7, \Di^\ 
markedly drops with increasing transit angle UTd when 
the latter exceeds unity. This effect is not important if 
the frequency-dependent GNR-FET transconductance is 
mainly limited by the electron delay under the top gate, 
i.e., if T ^ Td- However, at r ^ r^, the net transcon- 
ductance is determined by the product of \Gui\ and \Di^\ 
corresponding to curves shown in Fig. 4 and Fig. 7, re- 
spectively. 



VII. CONCLUSIONS 

(1) We proposed an analytical device model for a GNR- 
FET based on a heterostructure which consists of an ar- 
ray of nanoribbons clad between the highly conducting 
substrate (the back gate) and the top gate. 

(2) Using this model, we derived explicit analytical for- 
mulas for the GNR-FET transconductance as a function 
of the signal frequency in wide ranges of the collision 
frequency of electrons and the top gate length. 

(3) The variation of the GNR-FET high-frequency 
characteristics due to the transition from the ballistic and 
to strongly coUisional electron transport was traced. It 
was shown that the threshold frequency can nonmono- 
tonically depend on the ballsitic parameter. 



VIII. ACKNOWLEDGMENT 

The work was supported by CREST the Japan Science 
and Technology Agency, CREST, Japan. 



[1 

[2; 
[3; 
[4; 
[5; 
[6; 
[7: 

[lo; 



B. Obradovic, R. Kotlyar, F. Heinz, P. Matagne, T. Rak- 

shit, M. D. Giles, M. A. Stettler, and D. E. Nikonov, [11] 

Appl. Phys. Lett. 88, 142102 (2006). 

Z. Chen, Y.-M. Lin, M. J. Rooks, and P. Avouris, Physica [12] 

E 40, 228 (2007). [13] 

K. Wakabayshi, Y. Takane, and M. Sigrist, Phys. Rev. [14] 

Lett. 99, 036601 (2007). 

Y. Quyang, Y. Yoon, J. K. Fodor, J. Guo, Appl. Phys. [15] 

Lett. 89, 203107 (2006). 

G. Liang, N. Neophytou, D. E. Nikonov, and M. S. Lund- [16] 

Strom, IEEE Tran. Electron Devices 54, 677 (2007). 

G. Fiori and G. lannaccone, IEEE Electron Device Lett. [17] 

28, 760 (2007). 

V. Ryzhii, M. Ryzhu, A. Satou, and T. Otsuji, J. Appl. [18] 

Phys. 103, 094510 (2008). 

S. Briggs and J. P. Leburton, Phys. Rev. B 38, 8163 [19] 

(1988). 

S. Briggs, B. A. Mason, and J. P. Leburton, Phys. Rec. [20] 

B 40, 12001 (1989). 

G. Pennington and N. Goldsman, Phys. Rev. B 68 



045426 (2003) 

J. Guo and M. Lundstrom, Appl. Phys. Lett. 86, 133103 

(2005) 

W. Shocklcy, J. Appl. Phys. 9, 635 (1938) 

S. Ramo, Proc. IRE 27, 584 (1939). 

V. Ryzhii and G. Khrcnov, IEEE Trans. Electron Devices 

42, 166 (1995). 

V. Ryzhu, A. Satou, M. Ryzhii, T. Otsuji, and 

M. S. Shur, J. Phys. Cond. Mat. (in press). 

V. I. Ryzhii, V. A. Fedirko, and I. I. Khmyrova, Sov. 

Phys. Scmicond. 18, 527 (1984). 

V. I. Ryzhu, N. A. Bannov, and V. A. Fedirko, Sov. Phys. 

Semicond. 18, 481 (1984). 

V. I. Ryzhii and V. A. Fedirko, Sov. Phys. Semicond. 18, 

691 (1984). 

M. Shur, Physics of Semiconductor Devices (Prentice 

HaU, New Jersey, 1990). 

S. M. Sze, Physics of Semiconductor Devices (Wiley, New 

York, 1981). 



